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Abstract 

A family of orbiting resonances in molecular scattering is globally described by using a 
single pole moving in the complex angular momentum plane. The extrapolation of this pole 
at negative energies gives the location of the bound states. Then a single pole trajectory, that 
connects a rotational band of bound states and orbiting resonances, is obtained. These complex 
angular momentum singularities are derived through a geometrical theory of the orbiting. The 
downward crossing of the phase-shifts through 7r/2, due to the repulsive region of the molecular 
potential, is estimated by using a simple hard-core model. Some remarks about the difference 
between diffracted rays and orbiting are also given. 

1 Introduction 

After a long search, about two decades ago it became possible to observe orbiting resonances (or 
quasi-bound states) in molecular beam scattering [|^, ^, |3|, ^, ^. Although these phenomena can 
be quite naturally interpreted in the framework of scattering theory , their semiclassical nature 
calls for a more refined analysis, mainly geometrical, which has been partially performed by several 
authors, and notably by Berry |0, ^, |9|. In this spirit, it has been also advocated the use of the 



complex angular momentum plane polology, especially by Bosanac [10, 11, |12|], Connor [ll3| , [IJ] 
and Nussenzveig (see in particular his book |jl^ and the references quoted therein). In the latter 
approach, however, some problems remain, and a detailed phenomenological analysis is in part still 
missing. 

In the standard complex angular momentum approach one generally considers a given class 
of potentials, and then explores the analytical properties of the partial-waves along with their 
asymptotic behaviour in the complex angular momentum plane. If these analytical properties and 
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asymptotic behaviours allow for using the Watson transformation [^], then a resummation of the 
partial wave expansion can be performed, and the total scattering amplitude can be represented 
in terms of poles along with a background integral. The main limit of this method consists in 
the fact that it works only for a rather restricted class of potentials, notably for the Yukawian 



class |17], which is not very important in molecular scattering. Nussenzveig |18, |T^, in a 
series of very significant papers, extended the complex angular momentum method with particular 
attention to optical problems. In this context, the hard-core potential has been treated in detail, 
and, in addition to the Watson resummation, also the Poisson transformation has been used (see 
also Ref. [^]). The latter seems to be more promising for semiclassical scattering theories like the 
one that is of interest in the present paper. However, it must be noted that molecular scattering 



is usually described in terms of Lennard-Jones or Morse-type potential |21]; it is worth reminding 
that singular potentials of the type g^r~'^, n > 2, present an infinite number of poles in the first 
quadrant of the complex angular momentum plane, and, in such a situation, one is forced to use a 
modified Watson transformation. As far as we know, this approach has been rigorously proved in 
the case g'^r~^ by Dombey and Jones [^], and then conjectured for g'^r~^, n > 4. Furthermore, 
several authors, and notably those from the Uppsala school |2^, |2j, |2^, have applied the phase- 
integral method to derive locations and residues of the poles in the complex angular momentum 
plane, for singular potentials of the form V{r) = g'^r~"'f{r), (n > 4). It should be noted that the 
primary goal of all these researches consists in substituting to the partial-wave expansion, in a 
region where it is slowly convergent, more rapidly convergent representations. The latter then split 
the scattering amplitude into a sum of pole contributions and one or more integral terms which are 
usually evaluated by means of the saddle point method. 

Our viewpoint is quite different. We think that also in the region where only one partial-wave 
is dominant a complex angular momentum representation can be useful if we want to interpolate 
the various orbiting resonances with a single moving pole, whose extrapolation at negative energies 
can then give the bound states. With this in mind, we focus our attention on the poles of the 
scattering amplitude in the complex angular momentum plane, and we neglect, to a certain extent, 
the background integral, as will be explained below. To this purpose, a geometrical method of 
derivation of the complex angular momentum poles is implemented: instead of studying the ana- 
lytical and asymptotic properties of the partial-waves for a specific class of potentials, we analyze 
the particle motion (and specifically the orbiting) in a Riemannian space whose metric is induced 
by the potential. 

In view of the fact that we are treating a semiclassical scattering problem, the concept of particle 
path (or trajectory) still conserves its meaning. We are therefore allowed to study the geodesies 
(i.e., the particle paths) in a space whose Riemannian metric is given by dx' dx^ , r? = 1-V/E, 
E being the energy and V the potential. The coordinates require a more detailed comment: in 
fact, one generally claims to solve the dynamical problem in the whole space by using a unique 
coordinate system for all the points of the space being considered. Here we prefer to introduce 
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local coordinates that are appropriate to describe locally the particle trajectories, and then to 
use different coordinate systems elsewhere in space. We are thus led to consider the Jacobian 
of the transformation connecting the local coordinates to the ambient space coordinates. The 
transformation is singular in those domains where the Jacobian vanishes. We can then work out 
the problem by means of two geometrical tools: geodesies, which describe locally the particle paths, 
and Jacobians. The global solution emerges by patching up geodesic segments which represent local 
solutions. We will show in section || that by the sole use of these geometrical ingredients is possible 
to obtain a representation of the scattering amplitude in terms of complex angular momentum 
plane singularities (poles). 

However, this representation is not complete, and can only be regarded as an approximation; 
in what follows we shall discuss pro and con of this representation. First of all we obtain an an- 
gular distribution in terms of Legendre functions (not polynomials) of complex index, that present 
a logarithmic singularity forwards (i.e., at null scattering angle). This fact suggests that a com- 
pensating background term is missing in a representation which uses only poles. However, it is 
possible to project the scattering amplitude (represented only by poles) on the partial-waves, and 
the result can be regarded as a faithful representation, at least for angular momentum sufficiently 
small. Furthermore, the consequent integral (total) cross-section is finite, since the singularity of 
the Legendre functions involved is only logarithmic. On the other hand, the advantage of this 
representation is that the entire sequence of orbiting resonances can be fitted by using a single 
pole. The latter, while moving as a function of the energy, is able to describe in an ordered way 
the passage of the phase-shifts of various angular momenta through vr/2, that is the value of the 
phase-shifts corresponding to the resonances. Therefore, a global view and explanation of the 
orbiting resonances is obtained. 

Furthermore, these moving poles produce bound-states at those values of the energy where 
the angular momentum (which is real at negative energy) takes integer values. We thus obtain a 
trajectory of the complex angular momentum pole as a function of the energy, which gives bound- 
states at negative energy and resonances at positive energy. Obviously, the bound-states are not 
observed directly by scattering experiments. Their existence is, however, revealed by the fact that 
the phase-shift 5i at zero-energy is given by n/vr (as explained by the generalized Levinson theorem) 
where n/ is the number of bound-states with angular momentum equal to I (see refs. [p6| , [27[| , and 
sections ^ and These non-resonating phase-shifts decrease at higher energy passing through 
7r/2 in downward direction. This downward crossing of 7r/2 does not correspond to a resonance, 
but it is rather due to a repulsive potential. We can, however, fit this part of the phase-shifts by 
regarding them as the ones produced by a repulsive hard-core in a sense that will be explained in 
detail in section ^. 

The paper is organized as follows: in section |2| we outline a geometric theory of the orbiting, and 
obtain an approximate representation of the scattering amplitude in terms of moving poles (complex 
angular momentum poles) . In section ^ we derive a fitting formula for the phase-shifts that accounts 
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also for the role of the non-resonating phase-shifts. Section ^ is devoted to the phenomenological 
analysis. Finally, in section |5| we draw some conclusions, and briefly discuss the difference between 
orbiting and diffracted rays, and, accordingly, between resonances and diffraction. 



2 Outline of the geometrical theory of the orbiting 

Let us start from the Helmholtz's equation 

Au + k'^n'^u = 0, ^^2 = 1-^^, (1) 

where E is the energy and V the potential. Now, we look for a solution of equation (||) of the 
following form: 

n(x,fc) = y"^(x,/3)e''=*('''^)d/3, {^eR^). (2) 

The principal contribution to u{x,k), as k ^ oo, corresponds to the stationary points of ^, in 
the neighborhoods of which the exponential exp(ifc$) ceases to oscillate rapidly. These stationary 
points can be obtained from the equation d^{'x., (3)/d/3 = (assuming that 5^$(x, /9)/(9/3^ / 0). 
Then, the following asymptotic expansion of as — > oo, is valid: 

u(..t).e'»(-''o)f (3) 

m=0 ^ ^ 

where /3o is the unique stationary point of $ at fixed x. The leading term of expansion (^) reads 

n(x,A;) = ^o(x)ei'=*('''^°), (4) 

where 

-l/2\ 



^o(x) = A(x,/3o) 



exp < i-sgn — T ] } ■ (5) 



a/32 J 



/ /3=/3o 

The leading term (that hereafter will be written as Ae^^^ , omitting the subscript zero) can be 
substituted into equation (^) and then, by collecting the powers of (ifc) and nulling their coefficients, 
we obtain two equations: 

a) the eikonal equation 

mf=n' = l-^- (6) 

b) the transport equation 

V ■{A^V^) = Q. (7) 

The physical meaning of the transport equation (0) is that the probability current density is con- 
served. In order to exploit the geometrical content of the eikonal equation, it is convenient to 
replace the optical index v? with a Riemaniann metric tensor gij. Then, the eikonal equation can 
be rewritten as follows: 
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{x*} representing a system of local coordinates. Accordingly, the transport equation will be rewrit- 
ten in the following form: 

where g = |det {gij)\- Now, we consider the tube composed by the trajectories which describe only 
circular orbits and assuming, at this stage, that there is no leakage from the tube and, accordingly, 
there is no attenuation in the ray tube intensity. The angular coordinates used are 9 and (p and, for 
the sake of simplicity, we set the radius of the orbit to 1; then we have: gn = 1, 322 = sin^ 6, gi2 = 
g2i = 0. From Eqs. (^) and (|9|) we get 

(10) 

sme\A^^] \ =0, (e/nvr), (11) 





2 




) - 




1 


J" 


sin 9 





de J 

where we supposed that both the phase <I> and the amplitude A do not depend on (j). 

From equation ( p!o[ ) we have: $ = ±9 + const. . Next, substituting equation ( [lO| ) into equation 
(|lll) we finally obtain (see formula (^)) 

u{9, k) = -^^= e^'"' , (M nTT, n = 0, 1,2,...), (12) 
V I smt^l 

where the terms e^^^^ represent waves traveling in counterclockwise sense [e'^^^) or in clockwise sense 
^^-ik9y j^^ ^Yiis point let us note that approximation ( |l2|) fails at 9 = rni {n = 0,1,2,...): i.e., 
this approximation is not uniform. In particular, we are obliged to find the connection formulae 
relating the traveling waves when they cross the antipodal points 9 = 0, and 9 = tt. These 
connection formulae can be established by the use of the Maslov indexes. Here we simply give the 
results (the interested reader is referred to refs. |28, p9|]): 



i) for the counterclockwise orbiting wave we have, after the crossing through the point 9 = 0, 

after the 
crossing 

e"'' . e'^'^^-i) ; (13) 

ii) for the clockwise orbiting wave we have, after the crossing through the point 9 = 0, 

after the 
crossing 

e-'^' . e-^^'^^-i) . (14) 



Therefore, for each complete tour, both the counterclockwise and the clockwise oriented rays cross 
the antipodal points ^ = and 9 = n, and accordingly, in both cases the amplitude acquires a 
factor (—1). 
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We have seen that at the antipodal points the amphtude becomes infinite if we use the 9- 
representation. But in the Maslov construction use is made of the pe-representation (pQ being the 



variable conjugate to 6), and in the pg-representation the intensity of the ray-tube is finite |29]. 

In the case of unstable states the probability current density is not conserved and a tunnelling 
through the centrifugal barrier can, indeed, produce an attenuation of the intensity of the ray 
tube. The phenomenon can be represented by picturing trajectories which leave tangentially the 
circular orbit and emerge at large distances. In this process the whole space is involved and 
we must analyze the Jacobian J of the transformation relating the local coordinates, which are 
appropriate for describing the orbiting, to the coordinates of the ambient space, that is the whole 

space equipped with the euclidean metric. The leakage from the tube of circular trajectories 
can be described as follows. A beam of particles enter the interaction region, describe circular 
orbits of radius R, and, after a certain number of revolutions can emerge and leave tangentially the 
interaction region. In atomic scattering, the latter is not well-defined because of the long tail of the 
potentials. However, we can regard as the domain of interest for the orbiting process, the region 
delimited by the centrifugal barrier from one side, and by the hard-core (or more precisely by the 
high order singularity of the potential at the origin) from the other side. It is precisely in this 
domain that the orbiting resonances take place; the lifetime of these states is then related to the 
probability of tunnelling across the centrifugal barrier. The ambient space is properly described by 
the cartesian coordinates {x, y, z), the z-axis being directed along the incident beam and positively 
oriented in the direction of the outgoing particles. The local coordinates (^0)9^0)''") are defined as 
follows: 00 is the angle measured along the meridian circle from the point of incidence of the ray; 
(f)o is the azimuthal angle, and r is a parameter along the ray tangential to the circular orbit. We 
can then evaluate the Jacobian (see also Ref. |28||): 



J = = riRcosOo -rsmOo). (15) 

d(6'o,0o,T) 

The domain where the Jacobian vanishes is composed by: 
a) the surface r = 0; 

(3) the semi-axis (z > 0, a; = 0, y = 0) represented by r = r = i?cot ^o- 
When the Jacobian vanishes the eikonal approximation fails (for a detailed discussion of this very 



delicate mathematical point see Ref. |28]); furthermore, on the surface r = there is an attenuation 
of the intensity flux. We shall return below on this point. On the other hand, we can still use the 
transport equation written in the following form: 

iA(,^.) = „, (10) 

whenever J ^ 0, and the current density is conserved. In particular, in the asymptotic region we 
can write: A = const. 
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Now, since r = \Jr'^ — i?^, r being the radial variable (r^ = + + 2^), then r tends to r as 
r — > 00. Accordingly, we obtain: \fj —>■ ir^/sm 60 as r —>■ 00. At this point, in order to evaluate 
the contribution to the wavefunction, at large values of r, of the various classes of trajectories, we 
must distinguish between counterclockwise and clockwise oriented orbits. We start considering the 
contribution of the counterclockwise trajectories, that leave the circular orbit without completing 
one tour, and that cross the z-axis only once. We have, for large values of r, 

ik(R9++r) ikr ikRe+ 

{0 < 9q < tt). Formula ( |17D requires some comments. Hereafter we shall use a superscript (+) for 
all the quantities that refer to counterclockwise trajectories, and a superscript (— ) for all what refers 
to clockwise oriented rays; the subscript (0) in the notation of u refers to the fact that uq accounts 
only for the trajectories that have not completed one circular orbit. The factor a^k^O^) describes 
the leakage of the particles by the tunnelling effect. The phase $ (see formula (^)) now reads: 
^{6q ,r) = R9q + r ~ R9q + r, and the factor e~"^/^ can be regarded as the Maslov phase-shift 

r — >oc 

due to the crossing of the trajectory across the z-axis, where the Jacobian J vanishes (see formula 
(|l^)). We can now observe that the angle 9q coincides with the scattering angle 9s (see figure 
HI); we can then evaluate the contribution to the scattering amplitude due to the counterclockwise 
trajectories which have not completed one orbit. We obtain, with obvious notation: 

p-in/2 JkROs 

/+ (k, 9s) = aik, 9s) , (0 < 0. < vr). (18) 

^ ' V sm 9s 

Similarly, we can evaluate the contribution to the scattering amplitude of a clockwise trajectory 
which has not completed one tour. By observing that 9q = 2tt — 9s, we have 

ikR{2-K-es) 

f- {k,9s) = {-l)a{k,9s)^==^, {0 < 9s <7r). (19) 
^ ' V sm 9s 

Notice that the factor (—1) is due to the product of two Maslov phase factors, in view of the fact 
that the clockwise trajectory crossed the z-axis twice, even if it has not completed one orbit. 
Adding /^^^ to f^y we obtain the contribution of all the rays that have not completed one orbit: 

ikRBs _ ^JkR{2TT-es) 

fio){k,es) = /(+) + /(o) = -Mk,Os) ^— , {0<9s< vr). (20) 

We must now take into account the contribution of all those trajectories that describe more than 
one complete circular orbit before emerging. Let us assume that they do m [m E N) complete 
tours. The angles 9^^ are related to the scattering angle 9s as follows: 



O.m 



9s + 2TTm, (m = 0,1,2,...), (21) 



= 27r-0, + 27rm. (22) 
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We have, for < < vr, 



°o ^ikROs _ jgiA:i?(27r- 



f{k, 9,) = -i V i-ir e'^""'"' a{k, 9,) . (23) 

Once again the factor (—1), at each m, is due to the product of two Maslov phase-factors, corre- 
sponding to the fact that both the counterclockwise and the clockwise rays cross the z-axis twice 
for each complete circular orbit. The term a{k, 9s) can be regarded as a damping factor that gives 
the attenuation of the tube of circular orbits due to the tunnelling through the centrifugal barrier. 
At fixed k we can thus represent the leakage of the tube of orbiting particles by the use of an 
exponential of the following form: exp(— 7i?0^^). Accordingly, the damping factor a{k,9s) can be 
split in the product of two terms as follows: g{k) exp(— 7/?^^^). Coming back to the scattering 
angle 9s, we rewrite expression ( p3D as follows (see also formulae (^Tj) and (p^)): 

°o ifies _ i iM(27r-e,) 

f{k, 9s) = -ig{k) ^ (-1)- e'2--A< ^= , (24) 

(0 < < vr), where /i = R{k + 17). Next, we use the following expansion: 

-. 00 

= e'^^ V (-I)'" , (Im/u > 0). (25) 

2 cos TTU 

m=0 



By using equality (Eq) we can rewrite formula (p3) as follows: 



f{k, 9s) = -g{k) e-/^ , (26) 

2 cos vr^vsm 9s 



{0 < 9s < it). The r.h.s. of formula ( p6| ) contains the asymptotic behaviour of \J2'K[i times the 
Legendre functions P^_;^/2(~ cos 6*^), for 00 and |//|(vr — 9s) » 1 (see Ref. fS^). Then, 

writing i'^_i/2(— cos 9s) in place of its asymptotic behaviour, we have for \fi\ 00 and < < tt: 



2 cos TTfl 

Finally, by putting: /i — 1/2 = A, we obtain, for < 9s < it, 

f{K9s) = g{k) ^V2ATT ^^^:^°f-^ (28) 
2 smvrA 

It is then convenient to rewrite formula (^) in the following, more conventional, form: 

P\(F\{— cos 9) 

f{E, 9) = G{E) ' , {0<9<7t), (29) 

sm7rA(i/) 

where the scattering angle is now simply denoted by 9, and G{E) = ^y/^e''^/^g{k)y/2\ + 1 is the 
factor referred to the energy E. The r.h.s. of formula ( p9D is singular at 9 = 0, where P\[E)i~ cos 9) 
presents a logarithmic singularity [^]; therefore, representation (|2^ ) is only an approximation, since 
a background compensating term that makes the amplitude finite and regular at every value of 9 
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is neglected. Nevertheless this approximation is very useful because it can describe a sequence of 
resonances and bound states. In fact, if Re A is integer and ImA << 1, then the denominator 
sin7rA(i?) is close to zero and, accordingly, we observe a bump in the cross-section: i.e., we have an 
orbiting resonance. Since A is a function of the energy, representation ( |29| ) describes a sequence of 
resonances which are obtained whenever Re A crosses integer values, while Im A remains very small. 
As we shall see in section ImA tends to zero as E tends to zero. For negative energy ImA = 
and we have bound-states instead of resonances whenever Re A crosses integer values. At the end, 
we obtain a pole trajectory connecting bound-states to resonances as it will be shown in section ^. 



3 Orbiting Resonances and Bound— States 

In this paper we focus our attention only on those effects that contribute to the more relevant 
component of the spectrum: the orbiting resonances and the bound-states. However, they do not 
exhaust the entire phenomenology: in particular, there are also the "direct component" of the 
collision which is essentially composed by direct reflections, and the forward glory contribution. 
The former component does not require particular explanations; instead, it is worth spending a 
few words about the forward glory, especially in relation with the so-called "glory undulations". 
The latter effect is produced by those paths that, moving through attractive and repulsive regions 
of the potential, undergo equal positive and negative deflections and finally emerge at zero angle. 
Then, the interference between such paths and the forward-diffracted trajectories with large impact 
parameter give rise to forward glory oscillations that appear as undulations (glory undulations) in 
the total cross-section ITSl, |2l|l. The total number of oscillations can be related to the number of 



bound-states in the potential [21|. This relationship is one of the motivations of the great attention 
devoted to the "glory-undulations" by several authors ||3^. In our analysis we do not consider the 
"direct component" of the scattering, which can be globally described by a "background-term", 
which we neglect. 

More generally, let us briefly recall that, in classical scattering theory, if the deflection function 
goes through zero (or through a negative multiple of vr for a non-zero value of the impact parameter) 
the differential cross-section diverges like (sin©)"^ (B being the classical deflection angle). The 
divergence can occur either in the forward direction Q = —2mT (n = 0, 1,2, ... ; forward glory), 
or in the backward direction Q = — (2n + 1)tt (n = 0, 1, 2, . . . ; backward glory). Returning to the 
geometrical theory of the orbiting outlined in the previous section, we see that, up to formula (p6|), 
f{k,6s) presents singularities at 9s = and 9s = vr, which are due to the factor (smOg)^^^'^ ■ Next, 
passing from formula ( |2^ ) to representation ( p7[ ) we introduce the Legendre functions Px{—cos( 



which are logarithmically singular at 9s = 0, but regular at = vr. The smoothing of these 
singularities corresponds to the passage from the classical to the semi-classical approximation. As 
we have already noted several times, the obtained representation is far from being complete, as it 
does not represent the direct components of the collision: diffraction, direct reflection and forward 
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glory contribution |15|. 

Hereafter we will show that representation (29), in spite of the limits mentioned above, can, 
nevertheless, represent an ordered sequence of orbiting resonances. With this in mind, we project 
the amplitude ( p9[ ) on the Z-th partial wave, obtaining: 

^ e^-^Si - 1 ^ G{El 1 

2ik TT (a + i/?-0(a + i/? + / + l) ' 

where a + i/3 = A. Next, when the elastic unitarity condition may be applied, we get the following 
relationship among G, a and /?: 

G{E) = (3 {2a + 1) . (31) 

From ( pO| ) and (^) we obtain 

. . -1 (3{2a + 1) 

' {[(/-Q)2+/32][(/ + a + l)2+/32]}l/2 ■ ^ ) 

Formula ( [3^ ) represents a sequence of orbiting resonances: i.e., an ordered set of phase-shifts that 
cross 6i = Tr/2 with positive derivative. 

The bound-states are not directly revealed by scattering experiments. We know, however, from 
the generalized Levinson theorem that (5/(0) = n/vr, where n/ is the number of bound-states with 



angular momentum /, including those states which are not admitted by the Pauli principle |27]. 
Hereafter we shall consider, however, atom-atom scattering where the Pauli excluded states are not 
involved. Furthermore, in the cases considered in the next section we shall deal with a sequence of 
phase-shifts whose zero-energy value is exactly vr. This implies that there is a sequence of bound- 
states with increasing value of the angular momentum /. We shall see that these bound-states 
can be connected with continuity to the sequence of orbiting resonances, obtaining thus a band of 
bound-states and resonances (rotational band), generated by a unique pole moving in the complex 
angular momentum plane. At this point we are faced with the problem of finding a representation 
which is able to describe either those phase-shifts that do not resonate but decrease from their 
zero-energy value, that is tt, and cross tt/2 downward, and also that part of the phase-shifts that 
decrease after a resonance. Since these phase-shifts are due to the repulsive component of the 
potential, we can assume, in a very rough model, that they are produced by a hard-core potential. 
In the subsequent section we shall see that this model, in spite of its roughness, is in reasonable 
agreement with the phenomenological data. 

First, let us recall that the phase-shifts produced by an impenetrable sphere of radius ro are 



given by [g3| 

. , _i Ji+i/2ikro) 

Si = tan — — — , 33 

where J/+i/2 and Ni_^_i/2 denote the Bessel and Neumann functions, respectively. It is easy to check 
that when kr^ = {I + l)vr/2, then the phase-shift 6i, given by formula (p3|), equals {—tt/2); but the 
phenomenological data (see the next section) indicate that the radius ro of the hard-core, which 



E. De Micheli, G.A. Viano - Jan 11/02 



11 



models the repulsive component of the interaction, varies with /. We can then write kri instead 
of kro, and relate this expression to the angular momentum [a{a + 1)]"*^/^ through the equality 
kri = Bi[a{a + 1)]^^^, where Bi is a phenomenological parameter. We then glue together the term 
describing the orbiting to the hard-core phase-shifts as follows: 

. . -1 P{2a + 1) , 4 t -1 Ji+i/2iBiHa + ^)]'^'') 

where Ai, as well as Bi, depends on /, and both can be regarded as fitting parameters. As we 
shall show in the next section, formula (^) is able to reproduce an ordered sequence of orbiting 
resonances, and also gives in the correct order the downward crossing of the phase-shifts 5i through 
tt/2. The parameter Ai in formula ( p^ ) accounts for the strength of the repulsive action which 
models the singularity of the potential at the origin. This strength decreases as the angular mo- 
mentum increases, in agreement with the fact that the radius of the hard-core increases with /. 
Prom the phenomenological analysis the dependence of Ai and Bi on / turns out to be nearly linear 
(see figs, ^(c) and (d)). 

It is interesting to notice that a similar approach has been applied also in nuclear physics, in 
particular in the a-a and a-^^Ca elastic scattering (see refs. p^ , |35||). In these works the hard- 
core-type interaction has been used in place of the more realistic, but less tractable, exchange-forces 
which, ultimately, give rise to non-local potentials. However, in this case the approximation is not 
very faithful; in fact, while the fits of the rotational resonances are very satisfactory, for what con- 
cerns the downward passage of the phase-shifts, the theory yields only the general trend, but lacks 
for details. It is currently in progress an analysis of the complex angular momentum singularities 
associated with non-local potentials which aims at obtaining a more accurate representation of the 
downward behaviour of the nuclear phase-shifts. 



4 Phenomenological Analysis 

In this section a series of numerical fits, performed by using formula (|3^), will be presented. The 
experimental data being considered concern the H-Kr system, and are taken from Ref. in 
particular, we refer to figure 3 of that reference. In formula (^) the quantities Ai and Bi are 
considered as free fitting parameters, and for what concerns the functions a{E) and P{E) we adopt 
the following parametrizations: 

a{a + 1) = 2IE + oq , (35) 
P = dia{E)VE + d2VE , (36) 

where / = fiR^ is the moment of inertia, /i is the reduced mass, R the interparticle distance, and E 
is the center of mass energy; the parameters di and d2, that from the analysis turn out to be nearly 
independent of /, and the intercept ao are to be determined through the fits of the experimental 
phase-shifts (see the figure legends for numerical details). 
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In figure §(a) the fits of the phase-shifts for Z = 4, 5, 6, which display orbiting resonances, are 
shown. The upward crossings of 6i = tt/2, i.e. the energies of the resonances, are reported in figure 
|2|(b) in correspondence of the integer values of a = 4, 5, 6 (filled dots), along with their interpolation 
with the line a{a + l) as a function E. The good agreement between resonance locations {E > 0) 
and the line a{a + 1) strongly supports the correctness of the linear parametrization (|35D. Since 
ao > 0, we can extrapolate, for E < 0, the location of the bound-states corresponding to / = 0, 1, 2, 3 
(open circles). 

In figures ^(a) and (b) the fits of the phase-shifts for I = 0,1 and I = 2,3 are shown. The 
high energy behaviour is well reproduced, whereas the Levinson value at E = (i.e., 6i = vr), 
which is clearly satisfied by the experimental data, is not reached by the fitting curves. On the 
other hand, it is quite tempting to notice that if we extrapolate these theoretical phase-shifts to 
the unphysical region at negative energy (in spite of the fact that the phase-shifts have physical 
meaning exclusively at non-negative energy), they get close to the value 5i = it approximately at 
the energies Eq, Ei, E2 and £^3, that closely correspond to the location of the bound-states that 
have been recovered from figure ||(b), and whose angular momenta are precisely / = 0,1,2,3. In 
figures ^(c) and (d) the values of the parameter Ai and Bi that result from the fits of figures §(a) , 
|^(a) and ^(b) are displayed as functions of /. Ai presents a nearly linear decrease with /, while Bi 
increases almost linearly with /. 

Concluding, we illustrate the relationship between the function /? and the widths F of the 
orbiting resonances. These quantities are related to the tunnelling across the centrifugal barrier, 
which is now a complex quantity. We have 

from which it follows 

r = '^"^ — ^ , (/ = moment of inertia). (38) 

Prom (38) T can be easily evaluated, since a{E), I3{E), and / are recovered from the fits of the 
phase-shifts; the values of F are respectively: / = 4, F = 4.69- 10~^ meV; / = 5, F = 4.45- 10~^ meV; 
l = QV = 1.30 meV. 



5 Final Remarks 

In conclusion, we want to stress the difference between orbiting and diffracted rays. The diffracted 
rays can be explained as follows: when a ray grazes a boundary surface, it splits in two: one part 
keeps going as an ordinary ray, whereas the other part travels along the surface. At every point 
along its path this ray splits in two again: one part proceeds along the surface, and the other one 
leaves the surface along the tangent to the surface itself. This splitting can be mathematically 
explained as due to the non-uniqueness of the Cauchy problem for geodesies in a Riemannian 
manifold with boundary (see Ref. |^8|). The damping factors can then be related to the curvature 
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of the diffracting body ||2^. Conversely, for what concerns the orbiting case, the attenuation 
of the current of orbiting particles is produced by the leakage of the ray tube of the particles 
which describe circular orbits, and which is due to the tunnelling across the centrifugal barrier. It 
follows that the nature and the extent of the damping factors are drastically different in these two 
phenomena. We can say that, while the diffracted rays are produced by edge effects, this is not 
at all the case for the orbiting phenomena. Edge effects are indeed relevant in optics where the 
diffracting body (i.e., the obstacle) is compact, but in particle scattering they can be observed only 
in those processes described by short-range potentials, e.g. decreasing exponentially in r, which 
can simulate boundary effects. This does not happen in molecular scattering where the potentials 
that describe properly the phenomena are of Lennard-Jones type, which decrease, for large values 
of r, only as an inverse power. Thus, we can say that the near-forward diffractive scattering by a 
potential of this type, that is having a long-range tail, is not due to edge diffraction, but is rather 
produced by the contributions of very large angular momentum paths that undergo small deflection 
caused by the tail of the potential (see Ref. ||l5[| ). 

Finally, it is worth remarking that in the orbiting phenomena and in the case of diffracted rays, 
the attenuation can be explained merely as a damping of the amplitude, while the phase remains 
real: in both cases we can properly speak of real rays. However, one is tempted to explain the atten- 
uation of the intensity of the flux by introducing a complex phase, and, accordingly complex rays. 
We believe that this approach is not at all convenient for describing orbiting and diffraction phe- 
nomena, where the attenuation of the amplitude can be evaluated through the transport equation. 
Conversely, complex rays emerge, as appropriate mathematical solutions of the eikonal equation, 
for representing the shadow of the caustic. Therefore, they play a relevant role in total reflection, 
where they describe the exponentially damped penetration into the rarer medium associated to 
the surface waves traveling along the boundary. In molecular scattering they are again present in 
rainbow phenomena, and precisely in the shadow of the rainbow caustic. In fact, we can show that 
in the neighborhood of a caustic the nature of the eikonal equation differs from the standard one: 
its characteristics, which are real in the illuminated region, become complex conjugate beyond the 
caustic (i.e., in the shadow of the caustic), and correspond to the complex rays. 
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A X 




Figure 1: Geometrical model of orbiting: the ray below (solid line with open arrows) travels ^q" 
radians in counterclockwise sense; the ray above (dashed line with filled arrows) travels in clockwise 
sense and crosses the axial caustic (i.e., the z-axis) twice before emerging in the direction of the 
scattering angle 9s- The angle 6q coincides with the scattering angle Og- 
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Figure 2: H-Kr system, (a): Observed phase-shifts Q and corresponding fits (see formula (|3^) 
for the partial waves I = 4,5,6 versus the energy of the center of mass E. di = 0.13 (meV)^"*^/^; 
d2 = —0.48 (meV)~^/^ (see formula (|36[)). I = 4: observed (o), fitted (dotted line). / = 5: observed 
(+), fitted (dashed line). / = 6: observed (□), fitted (dash-dotted line), (b): a{a + 1) vs. E 
(see formula (^)). / = 5.1 (meV)~^; oq = 15.8. The filled dots indicate the energy Ei of the 
resonance for the partial waves / = 4, 5, 6. The open circles indicate the energy values Ei of the 
bound states for / = 0,1,2,3 extrapolated from the line a{a + 1): £^0 = —1-55 meV, Ei = —1.35 
meV, E2 = -0.96 meV, E^ = -0.37 meV. 
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Figure 3: H-Kr system. Observed phase-shifts |5[ and corresponding fits (see formula (^)) for 
the partial waves / = 0,1,2,3 versus the energy of the center of mass E. For what concerns the 
meaning of £"0, Ei, E2 and £"3 see the legend of figure ^. The values of the fits for E < (thinner 
lines) must be understood as explained in section (a): Case / = 0, 1. 1 = 0: observed (o), fitted 
(dash-dotted line). I = 1: observed (+), fitted (dashed line), (b) Case / = 2,3. 1 = 2: observed 
(o), fitted (dash-dotted line). / = 3: observed (+), fitted (dashed line), (c): Ai vs. /. (d): Bi vs. 
I. See formula (|34 



